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C/0 ■ ABSTRACT 

Of 

While there have been multiple observational programs aimed at detecting linear po- 
larization of optical radiation emitted by ultracool dwarfs, there has been comparatively 
less rigorous theoretical analysis of the problem. The general expectation has been that 
the atmospheres of those substellar-mass objects with condensate clouds would give 

y—t . rise to linear polarization due to scattering. Because of rotation-induced non-sphericity, 

there is expected to be incomplete cancellation of disk-integrated net polarization and 
thus a finite polarization. For cloudless objects, however, only molecular Rayleigh scat- 
tering will contribute to any net polarization and this limit has not been well studied. 
Hence in this paper we present a detailed multiple scattering analysis of the polariza- 
tion expected from those T-dwarfs whose spectra show absence of condensates. For 
this, we develop and solve the full radiative transfer equations for linearly polarized ra- 
diation. Only atomic and molecular Rayleigh scattering are considered to be the source 

• i-H , of polarization. We compute the local polarization at different angular directions in a 

plane-parallel atmospheres calculated for the range of effective temperatures of T dwarfs 
and then average over the whole surface of the object. The effects of gravity and limb 
darkening as well as rotation induced non-sphericity are included. It is found that the 
amount of polarization decreases with the increase in effective temperature. It is also 
found that significant polarization at any local point in the atmosphere arises only in 
the optical (B-band). However, the disk integrated polarization-even in the B-band-is 
negligible. Hence we conclude that, unlike the case for cloudy L dwarfs, polarization of 
cloudless T-dwarfs by atomic and molecular scattering may not be detectable. In the 
future we will extend this work to cloudy L- and T-dwarf atmospheres. 



Subject headings: stars: low mass, brown dwarfs — polarization, extinction — scattering 
— stars : atmospheres 
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INTRODUCTION 



Brown dwarfs have masses sufficient to ignite deuterium burning but insufficient to enter 
into the hydrogen burning main sequence. Instead, they radiate away the gravitational potential 
energy produced during their formation and contraction. As an individual brown dwarf cools it 
progressively passes through spectral types ranging from late M through the L and T sequences. 
Helpful reviews of brown dwarf science may be found in Chabrier & Baraffe (2000), Burrows et al. 
(2001), Kirkpatrick (2005), and Marley & Leggett (2008). 



to be formed in the atmosphere of relatively warm L dwarfs ( 


Lunine, Hubbard & Marlev 


1986; 


Tsuii et al. 


1996; 


Lodders, K. 


1999; 


Burrows, Marlev & Sharp 


2OO0l: Allard et al. 


2001 


), in the 



atmosphere of mid- type and later T dwarfs, grains condense well bellow t he photosphere and 



are not an important opac i ty source ( Ackcr man &; Marlev 



20001 : Burrows et al. Il200ll : iKirkpatrick 



2005 



2001 



Tsuii 



2002 



Chabrier Baraffe 



Marley &; Leggett 1120081 ). Synthetic spectra of model 



atmospheres lacking condensates fit well with the o bserved spectra of field T dwarfs later than about 



type T4 (with effective temperature T e g- ~ 1200 K; IStephens et al. I (120091 )) suggesting the absence 
of condensates in the visible atmosphere of these objects. At the same time, inclusion of condensate 
cloud of various species explains the spectra of L dwarfs (e.g., Cushing et al. 2008). 



Sengupta &; Krishan I (120011 ) first predicted that scattering of light by the condensates in the 
photopshere of L dwarfs would give rise to a significant amount of linear polarization. This po- 
larization would of course cancel out when integrated over a spherical surface of the star. How- 
ever, analysis of high resolution spectra of L dwarfs confirms that these object s are rapid rotator 



with the projected rotational velocity ranging from a s high as 90 to 10 kms 1 ( Basri et al 



200 



i 

Mohantv & Basri II2OO3I ; IZapatero Osorio et al. 1120061 ; iReiners fe Basri II2OO8I ; iBailer- Jones 112004 ). 
A rapid rotation around its axis induces departure from sphericity in the shape of the object. As a 
consequence, the net polarization integrated over the disk does not completely cancel out. Subse- 
quently, linear polari z ation of several L dwarfs has been detected by various group of astron omers 
( Menard et al. 1120021 ; IZapatero Osorio et al. 



observation were analyzed and explained by ISengupta I ([2003) ; ISengupta fc Kwok 



2005 



Goldman et al. 



2009 



Tata et al. 



(200 



20091 1. The 
13) who as- 



sumed a single scattering atmosphere wherein the effect of atoms and molecules was neglected. 
These authors suggested that scattering by atoms and molecules would give rise to negligible 
amount of polarization and hence polarization of T dwarfs in the optical should be inconsequential. 
However their single scattering approach left room for some uncertainty because they considered a 
purely scattering medium while a proper approach would incorporate both absorption and multiple 
scattering. 

In this paper we present the results of our detailed investigation on the polarization of T-dwarfs 
caused by multiple scattering of atoms and molecules in cloud-less atmospheres. In the next section 
we discuss the atmospheric models. In section 3, we present in detail the equations of transfer of 
polarized radiation and in section 4 we present the numerical methods of solving the transfer 
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equation of polarized radiation in a plane-parallel medium. The local polarization is integrated 
over the rotationally distorted surface. We briefly describe the method adopted to calculate the 
rotational oblateness and the methods of integrating the polarization over a rotationally distorted 
stellar disk in section 5. We discuss the results in section 6 followed by our conclusion in the last 
section. 



ATMOSPHERIC MODELS OF T-DWARFS 



The known field T dwarfs have effective temperatures in the range 1400 > T c g > 550 K. Clouds 
influence the observed spectra of the warmest dwarfs, but cloudless models generally reproduce 
the spectr a of most T dwarfs w ith T e g < 1200 K, although this temperature is likely gravity 
dependent ((Stephens et al. II2009I ). For the study re ported here we consider ed only cloudless models 



appropriate for spectral types later than about T3 ([Stephens et al. 1 12009 ) . 



For the polarization study we used four radiative-convective equilibrium models which give 
the run o f temperature as a functio n of pressure through the atmosphere from the collection em- 
ployed by lSaumon Sz Mar ley I (120081 ) in their study of brown dwarf evolution. The selected models 
i nclude no cloud opacity, although cond ensates are included in the chemical equilibrium calculation 
(jFreedman. Marlev and Lodders 1 120081 ). A more detailed summary of the model atmospheres is 
presented in lStephens et al. I (|2009l ). The models shown here all have log g = 5 and solar metallic- 
ity The effect of variation in these parameters will be explored in a future study. Here we choose 
models with T e g = 1400, 1200, 1000, and 800 K. As noted above, the first is somewhat warm for a 
typical cloudless field T dwarf but we include it as a comparison point for future studies of cloudy 
L dwarfs with this T e g- . 

The four model temperature-pressure profiles described above were then used as inputs for the 
polarization radiative transfer model. The same atmospheric opacities employed in the calculation 
of the radiative-convective equilibrium models (jFreedman. Marlev and Lodders 1 120081 ) were used 
in both calculations. We note that in the calculation of these temperature profiles we computed 
fluxes in 180 separate spectral bins. Within ea ch bin we computed the radiative transfer 8 times, 



once for each of eight k-coefficient ga uss points (jMarlev et al. 1120021) . The d etails of the underlying 



opacity calculation are presented in iFreedman. Marlev and Lodders I (120081 ) . For the polarization 
study we followed this same procedure and computed monochromatic intensities at centers of the 
same 180 flux bins using the same opacities. Final computed quantities are the gauss weighted sum 
of the eight independent radiative transfer calculations at each of the 180 wavelength points. We 
note that while the resulting spectrum has somewhat low spectral resolution, polarization changes 
slowly with wavelength and this resolution is adequate for our purposes. 
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3. EQUATIONS OF TRANSFER OF POLARIZED RADIATION 

3.1. Representation of polarized radiation by Stokes parameters 

In order to formulate the equations of transfer in a gaseous medium the most convenient 
representation of pola rized radiation is by a set of four parameters called Stokes parameters. 



Chandrasekhar I (|1960l ) first introduced the Stokes parameters in the equation of radiative transfer 
with a slight modification of Stoke's representation. In an elliptically polarized beam, the vibra- 
tions of the electric and the magnetic vectors in the plane transverse to the direction of propagation 
are such that the ratio of the amplitudes and difference in phases of the components in any two 
directions at right angles to each other are absolute constants. A regular vibration of this character 
can be represented by 

6 = ^ sm(ut - ei), & = 6° ] siaM - e r ) (1) 

where and £ r are the components of the vibration along two directions I and r at right angles to 
each other, u the circular frequency of the vibration, and ^ , ^ , Q and e r are constants. If the 
principal axes of the ellipse described by (&,£ r ) are in directions making angles x and X + \ 11 to 
the direction I, the equations representing the vibration take the simplified forms : 

£ x = f (°) C os j3 sin ut, £ x+ i T = £ (0) sin (5 cos ut (2) 

where (3 denotes an angle whose tangent is the ratio of the axes of the ellipse traced by the end 
point of the electric (or magnetic) vector, and the numerical values of it lies between and ^7r. 
The sign of (3 is positive or negative according as the polarization is right-handed or left-handed. 

In equation ([2]), denotes a quantity proportional to the mean amplitude of the electric 
vector and whose square is equal to the intensity of the beam: 



I=[^)f = [^f + [^ = I l + I r . (3) 

Following the representation given in equation ((2J) one obtains for the vibrations in the I and r 
directions 

£l = (cos (3 cos x sin ojt — sin (3 sin x cos ut) (4) 

and 

£ r = (cos (3 sin x sin ut — sin (3 cos x cos ut). (5) 
The intensities 1\ and I r in the directions I and r can be written as 

7, = [^ (0) ] 2 = 7(cos 2 cos 2 x + sin 2 (3 sin 2 X ) (6) 
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and 

I r = [£(°)] 2 = /(cos 2 (3 sin 2 X + sin 2 (3 cos 2 x). (7) 

From the previous equations it follows that if the regular vibrations representing an elliptically 
polarized beam can be expressed in the form given in equation dU), we can write the relations 

/ = [e; o) ] 2 + [^ 0) ] 2 = ^+^, (8) 

Q = [e{ 0) ] 2 - [^ 0) ] 2 = / cos 2/3 cos 2 X = h - I r , (9) 

U = 2^ (0) ^ 0) cos(e* — e r ) = I cos 2(3 sin 2 X = (Jj - I r ) tan 2 X (10) 



and 



V = 2f I (0) £<W sin(e/ — e r ) = I sin 2/3 = (J, - J r ) tan 2^ sec 2 X . (11) 

These are the Stokes parameters representing an elliptically polarized beam. It follows from equa- 
tion ([8j) to equation ([Tl]) that 

J 2 = Q 2 + u 2 + V 2 . 

Further, 

tan2 X = ^ 

and 

sin 2/3 



y /Q2 + [/2 + y2 

which give the plane of polarization and the ellipticity respectively. 

The nature of an arbitrary polarized light is completely determined by the intensities in two 
directions at right angles to each other (or equivalently, the total intensity / and Q = Ii — I r ) 
and the parameters U and V. The intensities I, Q, U and V are the general Stokes parameters 
representing light. 



3.2. Rayleigh scattering 



In order to incorporate Rayleigh's scattering into the radiative transfer equation. IChandrasekhar 



(|l960l ) modified it which can be stated as the vibrations representing the light scattered in a direc- 



tion making an angle with the direction of incidence is 



;(0) 



(|<7) 1/2 d 0) cos 9 sm{ojt - ei) (12) 
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and 



(o) 



;^)V^fsin(^- e2 ) 



(13) 



where the phase (ei, e 2 ) and the amplitude (£m°\£j°^), relations in the incident beam are maintained, 
unaltered, in the scattered beam. Here || and _L refer to directions in the transverse planes (of the 
incident and scattered light) parallel and perpendicular respectively to the plane of scattering. 
Accordingly, the parameters representing the scattered light are proportional to 



^[ej 0) ] 2 cos 2 e = ^/n cos 2 e, 



(14) 



3 rt(°)i2 3 t 



-a[2^ 0) ei_ 0) cos(ei - e 2 )] cos9 = -af/cosG 



and 



-a[2ej 0) ei 0) sin(ei - e 2 )] cos 9 = -aV cos 9. 



Therefore, denoting the incident light by the vector 

I = (I\\,I ± ,U,V) 
we can express the scattering intensity in the direction 9 by 

(a^)RIdw 



where 

/ cos 2 9 \ 

10 

cos 9 

\ cos 9 J 

R is called the phase matrix for Rayleigh scattering. 



(15) 
(16) 

(17) 



(18) 



3.3. The explicit form of the phase matrix for Rayleigh scattering 

In the formulation of the equation of radiative transfer, the radiation field at each point 
is characterized by the four intensities Ii(9,4>), I r (6,(f)), U{9,<p) and V(6,<f>) where 6 and (j> are 
the polar angles referred to an appropriately chosen coordinate system through the point under 
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consideration and I and r refer to the directions in the meridian plane and at right angles to it 
respectively. Therefore one writes 

I(9,<f>) = [Ii(0,<l>),Ir(9,<l>),U(0,<f>),V(9,<l>)]. 

The explicit form of the phase function for Rayleigh scattering in terms of 6 and 4> which is used 
in the transfer equation and which describes the angular distribution of the radiation field is given 
by (Chandrasekhar 1960): 



P(H,<I>; »',<//) 



Q[P<°W) + (1 - ^ 2 ) 1/2 (1 - ^P^W;^) 
+P( 2 W; M ',^)], 



(19) 



where 





/ 2(i-^)(iV) + /1 y 


M 2 





\ 


3 




1 








4 
























nn' ) 



( 4/x/i' cos(0' - <p) 2/zsin(0'-0) 







2/i' sin(<^' - (/)) cos(</>'-0) 

cos(0'-0) y 

/ fi 2 fi' 2 cos 2(4>' -4>) -jj 2 cos 2(0'- 



V 



p( 2 W; //,<//) 



~n' 2 cos 2(0' - 0) cos 2(0' - 0) 
/X// 2 sin 2(^' — 0) // sin 2(0' - 4>) 



\ 











// 2 // sin 2(0' 
-fi' sin 2(0' - 
//jii' cos 2(^' - 




\ 





and 



Q 



( 1 \ 
10 
2 
\ 2 J 

Here // = cos 6, fi' = cos 6*', and 0' denote the direction of the photon before and after scattering. 
Similarly for (j) and <p' . 

In case of the axial symmetry of the radiation field, it clearly requires that the plane of 
polarization be along the meridian plane (or, at right angle to it). Consequently, U = V = and 
the two intensities Ii and I r are sufficient to characterize the radiation field. 

With this consideration the phase matrix used to describe the angular distribution of photons 
that undergo Rayleigh scattering can be written as 



PGu,//) = 



2(l-/x 2 )(l- / x' 2 )+/xV 2 /j 2 



a 



(20) 
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3.4. Equations for the transfer of polarized radiation in plane- parallel geometry 

In the present work, plane parallel atmosphere has been considered. The transfer equation 
governing the intensities 1\ and I r in rest frame with the medium stratified into plane parallel can 
be written as 

4r ( ) ~ ( ) ( Z ) « - < l - <21> 

where \i = cos9(fi G [0, 1]) and 9 is the angle between the Stokes specific intensity vector and the 
axis of symmetry z; b(z) is the internal radiation source, Uq is the albedo for single scattering and 
r is the optical depth. 

In case of spherical geometry the equation of transfer for polarized radiation can be written as 

( Ii(ji) \ | l-/x 2 d ( 7,(x,/i,r) \ = 
dr y I r (n) J r dfi y I r (x,fj,,r) J 

{ ( Z ) - T L PM ( W> ) - (1 - "° )b(Z »} (22> 

where k is the mass absorption co-efficient for radiation of frequency v and p is the density of the 
material such that between the points s and s', 



t(s, S ')= r 

Js> 



upds. 



For an object having radius as small as that of Jupiter with geometrically thin atmosphere, 
the second term in the left hand side of the above equation which incorporates the curvature effect 
becomes negligible. Therefore, for solar type of stars and substellar-mass objects, the plane-parallel 
stratification is sufficient. 

The degree of linear polarization p is given by 

Ii(p) - I r ([i) 



4. NUMERICAL METHODS 



The work of Redheffer 



(1962) and iPreisendorfer I (|1965l ) on the interaction principle have 



been formalized by Grant &z Hunt ( 19691 ). Grant &z Hunt (1969a) with the i ntroduction of the 



i ntern al so urces which is cruci a l for t he stellar atmospheres. Following this work, I Grant fc Peraiah 



(|1972l ) and iPeraiah &; Grant I (|1973l ) developed a method to obtain direct solution of the transfer 
equations. This method is called the discrete space theory of Radiative Transfer. In the present 
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work the method due to lPeraiah fc Grant I ()1973l ) has been employed to solve the radiative transfer 
equations in their vector form for linear polarization. 

In this method, the entire medium is divided into N number of shells and it is assumed that 
the specific intensities l£ and I~, 1 are incident at the boundaries n and n + 1 respectively of the 
shell with optical thickness r. The symbols with signs + and - represent specific intensities of the 
rays travelling in opposite directions. If « represents the cosine of the angle made by a ray with 
the normal to the plane parallel layers in the direction in which the geometrical depth decreases. 
That is, 

I+[IM : < M < 1] 

and 

/"[/ n (-/i):0</x<l], 

J + represents the specific intensity of the ray travelling in the direction u and I~ represents the 
specific intensity of the ray travelling in the opposite direction. We select a finite set of values of 
fi(fJ-j ■ 1 < 3 < m; < m < fjL 2 < «3 ■ ■ ■ «m < 1) 

In{V>2) 



It 



\ Intern) ) 



and 



/ In(-»l) 
I~ = I n (-fMi) 

\ ■ In{ fJ"m) 

are m-dimensional vectors on Euclidean space. 

The incident intensity vectors are 1+ and I^>i- The emergent intensity vectors are I~ and I^+i- 
The emergent radiation field will have the contributions from the internal sources, say, S + (n + 1, n) 
and S~(n, n + 1) corresponding to the output intensity vectors I f J +1 and I~ respectively. 

We assume certain linear operators which reflect and transmit the incident radiation e.g., 
t(n,n + 1), t[n + l,n), r(n, n + 1) and r[n + l,n). These operators are calculated based on the 
physics of the medium. Then we can write the output intensities in terms of the transmitted and 
reflected input intensities together with the internal sources as 



In+i = t( n + !> n )In~ + r ( n , n + + S+(n + 1, n) 

I~ = r(n + l,n)/+ + t(n,n + + S _ (n,n + 1). 



(24) 
(25) 



The relationship given by equation (I24p and equation (1251) is called the Interaction Principle. Equa- 
tion (JMJ) and equation (|25|) can be written concisely as 



I 



n+l 

1,7 



S(n, n + l) 



It 



+ £(n,n + 1) 



(26) 
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where 



S(n, n + 1) 



t(n+l,n) r(n, n+1) 
r(n+l,n) t(n, n + 1) 



(27) 



If there is another shell wit h boundaries (n+1, n+ 2) adjacent to (n,n +1), interaction principle 
for this shell can be written as ([Grant fc Hunt 1 11969) 



4+2 I -= S(n + l,n + 2) ( V 1 )+E(n + l,n + 2) 



4+i 



(28) 



n+2 



where S(n + l,n + 2) is similarly defined as in equation (|27p . If we combine the two shells (n, n + 1) 
and (n + 1, n + 2) then the interaction principle for the combined shell is written as (for the thickness 
is arbitrarily defined): 



ln +_ 2 j = S(n, n + 2) I /n ] + E(n, n + 2). 

4 / \ 4+2 / 



(29) 



S(n,n + 2) is called the star product of the two S-matrices S(n,n + 1) and S(n + l,n + 2); and 
S(n, n + 2) can be written as 



S(n, n + 2) = S(n, n + 1) * S(n + 1, n + 2). 



(30) 



Equation (j29|) is obtained by eliminating 4~+i and I n+1 from equation ([26]) and equation (p8|) . We 
can write r and t operators for the composite cell as 



t{n + 2, n) = t(n + 2, n + 1)[J - r(n + 2, n + l)r(n, n + l)] _1 t(n + 1, n), 



(31) 



t(n,n + 2) = t(n,n+ 1)[J - r(n,n + l)r(n + 2,n + l)] _1 t(n + l,n + 2), 



(32) 



r(n + 2,n) = r(n + 1, n) + t(n, n + 1)[J — r(n + 2, n + l)r(n, n + 1)] 1 x 
r(n + 2, n + 1), 



(33) 



and 



r(n,n + 2) = r(n + l,n + 2) + i(n + 2, n + 1)[J - r(n, n + l)r(n + 2, n + 1)] 1 x 
r(n, n+1) 



E(n, n + 2) = A(n, n + 1; n + 2)E(n, n + 1) + A'(n; n + 1, n + 2)E(n + 1, n + 2) 
where I is the identity matrix and 



A(n, n + l;n + 2) 



t(n + 2,n + 1)[J - r(n + 2,n + l)r(n,n + 1)] _1 
i(n, n + l)r(n + 2, n + 1)[I - r(n + 2, n + l)r(n, n + 1)] _1 I 



(34) 



(35) 
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A'(n;n + l,n + 2) 

and 



J t(n + 2, n + l)r(n, n + 1)[I — r(n, n + l)r(n + 2, n + 1)] 1 
r(n, n + 1)[I — r(n, n + l)r(n + 2, n + 



S(n, n + 1) 



£+(rc + l,n + 2) 
£~(n, n + 1) 



Similarly, S(n + 1, n + 2) is denned. 



In order to obtain physical interpretation of the equations (|3ip - (|34|) we expand the operator 
inverse in a power series. For example, 



t(n + 2, n) = ^2 t k( n + 2, n), 



k=0 

t k (n + 2, n) = i(n + 2, n + l)[r(n, n + l)r(n + 2, n + l)] fe i(n + 1, n). 

This operator acts on intensities to the right and gives the contribution to /i_ 2 from J+. The term 
tk( n + 2, re) may be recognized as diffuse transmission from n to n+1, diffuse reflection from the 
layer (n, n + 1), k times in succession and finally diffuse transmission through (n + 1, n + 2). Thus 
i(n + 2, n) is the sum of contributions involving scattering of all orders k = 0,1, 2, ... , oo. A similar 
interpretation can be given for other operators. 

If we write S(a) (a to designate the cell) then 

S(a*/3) = S(a)*S(/3) (36) 

where a * (3 denotes the region obtained by putting the two cells a and f3 together. If the cells are 
homogeneous and plane parallel then 

a */3 = f3*a. (37) 

In general, star multiplication is non-commutative. However, star multiplication is associative. If 
we have to add several layers a, (3, 7, . . . then, 

[(a* (0*7) ★...)] = S[(a*/3)* 7 *...]. (38) 

If the medium is optically very thick then we can use what is known as 'doubling method'. For 
example, 

S(2 p d) = S(2 p - 1 d) * S(2 p - 1 d), (P = I, 2, 3, . . .) (39) 

which means that we can generate the S-matrix for a layer of thickness 2 p d in P cycles starting 
with S(d) rather than in 2 P cycles of adding the S(d)s one by one. For example if P = 10, then 
only a fraction 10/2 10 ~ 10 -2 of the computational work is needed to add 2 10 layers of thickness d. 

We expect the reflection and the transmission operators to be non- negative on the physical 
grounds that intensities are always non-negative. This condition will be satisfied only when the 
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optical thickness of the shell is less than certain value called the 'critical size' or r cr jj. If the optical 
thickness r of the shell in question is larger than t ct h then we can divide the shell into several 
sub-shells whose thickness r is less than T cr a and then use star algorithm to calculate combined 
response from the sub-shells whose total thickness is T. If, for example, we need the radiation field 
at internal points in the atmosphere, we shall have to divide the entire medium into as many shells 
as we need and calculate the radiation field at the N points in the medium. One can write down 
the interaction principle for each shell and solve the whole system of equations. 

The solution and (for any shell between shell 1 (at the outermost region) and shell N 
(at the innermost region) are obtained from the relations 

J+ +1 = r (l, n + 1)I- +1 + V+ +1/2 (40) 

and 

I- =t(n,n+l)I- +1 + V- +1/2 (41) 

with the boundary conditions = I (a)- The quantities r(l, ra+1), V^ +l and V~ +l are calculated 
by employing the initial conditions r(l, 1) = and Vy 2 = I + (b). The computation is done by the 
following recursive relations : 

r(l,n + 1) = r(n, n + 1) +t(n+ l,ri)r(l,n)[I — r(n + 1, n)r(l, n)]~ 1 t(n, n + 1), 

V n+l/2 = f ( n + ^ n ) V n-l/2 + S+ ( n + n ) + R n+l/2^' (™, " + 1), 

V n+l/2 = r(n+l,n)V+_ 1/2 +T n+1/2 Z-(n,n+l), (42) 

where 

t(n + 1, n) = t(n + 1, n)[I — r(l, n)r{n + 1, n)]^ 1 , 
r{n + 1, n) = r(n + 1, n)[I — r(l, n)r(n + 1, n)]~ , 
Rn+i/2 = t{n+l,n)r(l,n), 
T n +i/2 = [I ~r(n+ l,n)r(l,n)] _1 , 

t(n,n + l) = T n+1/2 t(n,n+l). (43) 

To calculate the radiation field at the internal points we proceed as follows : 

(1) Divide the medium into a number of shells (say N) with N+l boundaries as mentioned earlier. 

(2) Start calculating the two pairs of reflection and transmission operators r(n + 1, n), r(n, n + 1), 
t(n + 1, n) and t(n, n + 1) in each shell. If the optical thickness of any shell is larger than T cr a then 
apply star algorithm to use doubling procedure if the medium is homogeneous. 

(3) With the boundary condition that r(l, 1) = and V^ 2 = I + (b) and the r and t operators 
mentioned in (2) compute recursively r(l,n + 1), Vy 2 and t(n,n + 1) given in equation (f4lj) to 
equation (|43p from shell 1 to shell N i.e., from b to a. 
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(4) Next sweep back from a to b calculating the radiation field given in equation (|40p with the 
boundary condition = I~(a). 

We use a eight point Gaussian quadrature formula for the numerical integration over the 
angular points and hence /i and // are the zeros of the Legendre polynomials of degree 8 in the 
interval -1 to +1. The phase matrix is normalized to 1. 



5. ROTATION-INDUCED OBLATENESS AND DISK INTEGRATED 

POLARIZATION 



Once and I r {l^) are calculated in any local points, the polarization integrated over the 

surface of the object can be calculated which should be the observable quantity for point source 
of light. However, the net polarization would cancel out to zero because of symmetry if the the 
apparent disk of the object is perfectly spherical. Rapid rotation of the object makes it non-spherical 
and hence non-zero polarization arises when integrated over the surface. 



The oblateness of a rotating object has been discussed bv lChandrasekhar I (|1933l ) in the context 
of polytropic gas configuration under hydrostatic equilibrium. For a slow rotator, the relationship 
for the oblateness / of a stable polytropic gas configuration under hydrostatic equilibrium is given 
by 

n 2 m 



J R e 3 GM 



(44) 



where M is the total mass, R e is the equatorial radius, R p is the polar radius and O is the angular 
velocity of the object. C is a constant whose value depends on the polytropic index. For the 
polytropic index n = 0, the density is uniform and C = 1.875. This configuration is known as 
the Mac laurin spheroid. For a polytropic index of n = 1.0, C = 1.1399, which is appropriate for 
Jupiter (IHubbard Ill984l ). For non-relativistic completely degenerate gas, n = 1.5 and C = 0.9669. 
Comparisons with detailed structure models (D. Saumon, private communication) show that brown 
dwarf interiors can be adequately approximated by polytropes with 1 < n < 1.3 with the larger 
n being appropriate for higher gravities. For n=l and R e = Rj U p the oblateness varies from 0.001 
to 0.086 when the projected rotational velocity varies from V=10 to 90 kms -1 with the surface 
gravity g = 10 5 cms -2 and it varies from 0.00035 to 0.0287 with g = 3 x 10 5 cms -2 . Similarly, for 
n=1.5 and R e = i?j up , the oblateness varies from 0.0009 to 0.073 with g = 10 5 cms and it varies 
from 0.0003 to 0.024 with g = 3 x 10 5 cms -2 for the same variation in the rotational velocity. It 
should be noted in this context that at 1 bar pres sure level, the obla t eness / of Jupiter, Saturn 
and Earth are 0.065, 0.098 and 0.003 respectively. iBarnes &: Fortney I ([2003) used Darwin-Radau 
relationship and estimated the oblateness of the exoplanet HD209458b to be about 0.00285 whereas 
the polytropic approximation with n=l yields a value of 0.00296. 



Therefore, in the present work, we employ the polytropic approximation of Chandrasekhar with 
the polytropic index n = 1, to find out the rotation-induced oblateness because this is appropriate 



-14- 



for low mass dwarfs and provides an upper limit to the oblateness and hence an upper limit on the 
degree of polarization for a given rotational velocity. As stated before, the surface gravity of all the 
models is fixed at g = 10 5 cms -2 . 

We assume an axisymmetric shape of the object so that it has a rotational invariance around 
some axis. The radius of any point on the surface is given by 



R e 



[1 + {A 2 - l)/! 2 ] 1 / 2 



(45) 



where A is the ratio of equatorial radiu s to the polar r adius, i.e., A = 1/(1 — /). Using the spherical 
harmonic expansion around 6 and (f) (j Jackson I Il975l ). we find the total flux integrated over the 
surface as 



I 00 i-i 
F I = -Y^(2l + l)Pl{co B i) 

1=0 



[1 + (A 2 - 1)^ 2 ]V2 J_ x 



(46) 



and the polarized flux as 



oo „! 

F Q = 2vrVa 2 (/, 2)P, 2 (cosi) / 
1=2 J- 1 



Pi(n)d/j, 



[1 + (A 2 - l)/i 2 ]V2 J_ x 



[I^-Ir^flP^dfi 



(47) 



where 



a(l, m) 



(2/ + l)(/-m)! 



4ir(l + m) 



1/2 



(48) 



Pi{[i) is the Legendre polynomial and P™^) is the associated Legendre polynomial. The disk 
integrated linear polarization is p = Fq/Fj. 



The disk integrated polarization of a star is investigated by lHarrington Collins II I (| 19681 ) 
who assumed a Roche model wherei n the entire gravitationa l mass is considered to be at the 
center of the star. The formalism of lHarrington &: Collins II I (119681 ) although is valid for stars, 
is unlikely to be applicable for brown dwarfs and giant gaseou s planets because of the i r very 
different mass distributions. Nevertheless, we use the formalism of lHarrington &: Collins II I (J1968) 
and calculate the disk integrated polarization in order to obtain an order of magnitude check in 
the results with spherica l harm onic expansion method. However, the shape function x(£l,8) in 
Harrington &; Collins II I (|1968l ) is replaced by 



x(0, n) 



1 



[1 + (1M 2 -1)^ 2 ]V2 



(49) 



so that the shape of the object takes the form of an ellipsoid, as before, and the formalism becomes 
comparable with the spherical harmonic expansion method. In both the methods we consider 
R e = Rjup and the mass M is derived from the surface gravity which is fixed at 10 5 cms -2 . 
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6. RESULTS AND DISCUSSION 

The effective temperature of T dwarfs varies from about 550K to 1200K. Above 1200K, conden- 
sates of various species appears in the photosphere and transition from L dwarfs to T dwarfs takes 
place at about 1300K to 1400K. On the other hand, T dwarf with effective temperature less than 
800K are likely too faint for image polarimetry. Therefore, in the present work we have considered 
an effective temperature ranging from 800K to 1200K. However, for a comparison with L dwarfs of 
effective temperature 1400K which have condensates in their photosphere and hence should show 
detectable amount of polarization in the optical because of dust scattering, we have included the 
case for T e g = 1400-fT without condensates. This case provides a point of comparison to illus- 
trate how polarization differs with or without dust sc attering. The surface grav ity of typical field 



brown dwarfs varies from 10 5 to 3 x 10 5 cms 2 (e.g., ISaumon &: Marlev I ( 200S)). The departure 



from rotation-induced sphericity decreases with the increase in surface gravity (jSengupta k, Kwok 



20051 ). In order to maximize the amount of polarization and hence to obtain an upper limit on the 



degree of polarization of T dwarfs, we have fixed the surface gravity at g = 10 5 cms 2 . 

The numerical method used to calculate the local polarization provides unconditional stability 
for plane-parallel stratification of the atmosphere subject to the condition that each shell has less 
than the critical optical depth. This is achieved by dividing the shells with high optical depth 
into sub-shells with optical depth less than the critical optical depth. Usually, r < 0.1 serves 
the purpose. A convenient test of the efficacy of the numerical method is to study the case of 
conservative scattering. In a purely scattering medium, the physical system must neither create 
nor destroy energy. For this purpose we apply unpolarized incident radiation at the boundary of 
the inner radius of the shell and no radiation is incident at the outer shell. For a conservatively 
scattering medium, the total flux that is introduced at the inner boundary is found to be equal 
within nine-th decimal figure to the sum of the flux that comes out of the outer boundary and 
the backscattered flux at the inner boundary. The polarization is zero when the phase function 
P(/x, //) = 1 i.e., when isotropic scattering is considered. 

Finally, the total emergent flux given by 

F = F l + F T = J [/,(//) + I r (p)]ndfjf (50) 



matches exactly with the emergent flux obtained by iMarlev I ((2000) and IStephens et al. I ((2009) 
for scalar radiative transfer case. In figure 1, we present the emergent flux for condensate- free 
atmospheres with effective temperature ranging from 800K to 1400K. As it is well known, the 
emergent flux in the optical falls rapidly with the decrease in wavelength owing to the presence of 
alkaline elements (Burrows, Mar ley & Sharp 2000; Liebert et al. 2000). Therefore, T dwarfs are 
extremely faint at shorter wavelengths. 

Figure 2(a-d) shows the polarization profiles at different directions for different effective tem- 
peratures. The polarization is maximum when 6, the angle between the direction of the radiation 
field and the symmetry axis is 90°. However, the degree of polarization falls steeply with the 
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decrease in 9. At 9 = 90°, significant polarization is found even up to a wavelength of 1.5 //m 
for the coolest T dwarf (T e g = 800-KT). But at a slightly smaller value of 9, polarization becomes 
insignificant at wavelength longer than even 0.6 /im. The angular dependency of the polarization 
for different effective temperature is presented in figure 3. As the angle between the direction of 
the radiation field and the symmetry axis decreases, the anisotropy in the radiation field reduces 
drastically because of the fact that at smaller angles, mostly the unscattered photons emerged out 
and that the probability of scattering increases as 9 increases. The probability of scattering is 
maximum at 9 = 90°. 

Figure 2(a-d) also shows that the degree of polarization decreases as the effective temperature 
of the object increases. In fact, this is true, not only at the outermost boundary but also at any 
depth inside the atmosphere. In figure 4 we present the atmospheric depth dependence of the 
polarization at a particular wavelength (A = 0.6/xm) and at a direction \x = 0.02. The degree of 
polarization for T e g = 800-RT remains almost constant up to 1 bar of pressure level and then falls 
rapidly to zero. As the effective temperature increases, the degree of polarization falls rapidly to 
zero at a higher altitude or at a lower pressure level. The physical explanation of this feature 
becomes clear from figure 5 wherein the variation of the single scattering albedo ujq with respect 
to the atmospheric pressure P and temperature T is presented. For an object with T e g- = 1200i^, 
the scattering albedo is almost constant and its value is about one above a pressure level of 10 _1 
bar. Bellow this depth, it asymptotically falls to zero, uq becomes zero at deeper region as one 
goes from hotter to cooler objects. This means, the scattering albedo remains non zero up to much 
deeper region as the objects become cooler and hence contribution to the polarization originates 
from deeper region in cooler objects. As a consequence, the degree of polarization is higher in 
cooler objects as shown in figure 4. 

The polarization profile discussed above is calculated at a local point of the surface of the 
atmosphere stratified into plane-parallel geometry. This can be observable only if the object can be 
spatially resolved. However, a distant stellar object cannot be spatially resolved and it appears to 
be a point source of light to an earth based observer. Therefore only the disk averaged polarization 
is the observable quantity. However, if the apparent disk of a stellar object is perfectly spherical, 
because of symmetry, the net polarization averaged over the spherical disk would cancelled out. 



High resolution spectroscopic analysis of T dwarfs bv lZapatero Osorio et al. I (|2006l ); ldel Burgo et al 



(J2009J) shows that just like L dwarfs, these objects are also fast rotators. As a consequence, net 
non-zero polarization should arise when integrated over the apparent disk because rotation induces 
distortion in the stellar disk. In figure 6 and figure 7, we present the disk integrated polarization 
of T dwarfs at an edge on view, i.e., at an inclination of 90° with rotational velocity V=90 and 60 
kms -1 respectively. The degree of polarization is maximum when the inclination angle is 90°, i.e., 
at an edge on view to a n observer and it decrease s with the decrease in the inclination angle as can 



be seen from figure 8. lHarrington &: Collins II I (|1968l ) also reported the same feature. Figure 6 



and Figure 7 shows that the degree of polarization of T dwarf with any spectral type is non-zero 
only at wavelengths shorter than 0.6 //m. As the rotational velocity increases, the oblateness of 
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an object with the same mass, radius and polytropic index increases. As a consequence with in- 
creasing rotational velocity there is less cancellation of polarization over the surface, yielding higher 
polarization. At the same time, the effect of limb darkening and gravity darkening too increases 
and consequently the degree of polarization increases further. However, because of the fact that 
T dwarfs are extremely faint at wavelengths shorter than 0.6 /im, it would be extremely difficult 
to detect such a low polarization of T dwarfs, especially if the inclination angle is off the edge on 
view. 

The synthetic spectra calculated by using one dimensional plane- parallel atmospheric strat- 
ification match well with the observed spectra. However, a three dimensional approach is worth 
investigating in order to achieve better understanding of the effect of limb darkenning and the 
polarization integrated over the disk of a non-spherical object. 



CONCLUSIONS 



In the present paper, we have discussed the results of our investigation on multiple scattering 
polarization of cloudless T dwarfs. W e predict that in the ab sence of condensates in the photosphere 
of T dwarfs later than about type T3 ([Stephens et al. 1120091 ). non-zero polarization would arise only 
at wavelengths shorter than 0.6/iin- a spectral region where these dwarfs are exceptionally dark. In 
order to calculate the maximum possible polarization of T dwarfs, we have considered a polytropic 
equation of state with index n = 1 which is the minimum permissible value of n for representing the 
matter distribution in brown dwarfs. Consequently, the rotation- induced oblateness is maximized 
for a given rotational velocity. 

We find that even with a rotational velocity as high as 90 kms -1 , the maximum amount of 
disk integrated polarization that arises at an inclination angle of 90°, is of the order 10~ 4 , too 
small to be detected in the B-band. The degree of polarization would be greatest for the coolest 
ob jects and it would decrease with increasing e ffective temperature. The results support the claim 
bv lSengupta I (|2003l ); ISengupta &; Kwok I (|2005l ) that the polarization detected in the optical region 
(R- and I-band) of L dwarfs has negligible contribution from atomic and molecular scattering. 
Therefore, if linear polarization in R- or I-band of T dwarfs is detected in future, it will imply the 
presence of condensates in the photosphere of T dwarfs. Although a magnetic field in principle 
might produce polarization, T dwarf atmospheres are highly neutral (Gelino et al. 2002) and 
thus th e presence of magneti c field would not give rise to detectable amount of polarization in the 



optical ([Menard et al. 1 12002). Finally, polarization could play a crucial role in measuring the loss 



of condensates from brown dwarf atmospheres as we expect L- and early T-type dwarfs to show 
linear polarization in R- and I-band by dust scattering in the phototosphere whereas mid-type and 
later T dwarfs would not show any polarization at wavelength greater than 0.6/um. The detailed 
investigation of the polarization of L dwarfs by multiple dust scattering will be presented in a 
forthcoming paper. 
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Fig. 1. — The emergent flux as a function of wavelength for T dwarfs with different effective temper- 
ature. From top to bottom, the solid lines represent the flux for T e// =1400,1300, 1200,1100,1000,900 
and 800K respectively. 
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Fig. 2a. — Percentage degree of polarization along fi = 0.02 for T dwarfs with effective temperature 
800K, 1000K, 1200K and 1400K. 
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Fig. 2b. — Same as figure 2a but along fi = 0.4. 
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Fig. 2c— 



Same as figure 2a but along ^ = 0.76. 
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Fig. 2d.— 



Same as figure 2a but along ^ = 0.98. 



-26- 




Fig. 3. — Percentage degree of polarization at different fi = cos where is the angle between the 
direction of the radiation and the symmetry axis. The polarization is calculated at a wavelength 
A = 0.6 fim. 
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Fig. 4. — Percentage degree of polarization along fi = 0.02 at different atmospheric pressure level. 
The polarization is calculated at a wavelength A = 0.6 /im. 
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Fig. 5. — Atmospheric depth dependence of single scattering albedo u;q at a wavelength A = 0.6 
jum. The thick solid lines represent loq vs. atmospheric pressure P scaled in the left hand side axis. 
The dashes over the lines indicate the location of the photosphere. The thin solid lines represent 
the variation of loq with respect to the atmospheric temperature T scaled in the right hand side 
axis. 
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Fig. 6. — The disk integrated polarization at an inclination angle i = 90° and for a rotational 
velocity V = 90 kms -1 . The solid lines represent the percentage polarization calculated by using 
the spherical harmonic expansion method while the broken lines represent that calculated by using 
the modified Harrington-Collins method. From top to bottom the lines represent the disk integrated 
polarization of T dwarfs with effective temperature 800K, 1000K and 1200K respectively. 
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Fig. 8. — Disk integrated polarization at different inclination angles i. The polarization profiles 
shown here are calculated by using the spherical harmonic expansion method. 



